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Abstract 

Chemical and other species reaction theories involving thermodynamical equi- 
librium states characterized by a temperature parameter invariably utilize 
statistical mechanical equilibrium density distributions. Here, a definition of 
heat-work transformation termed thermo mechanical coherence is first made, 
and it is conjectured that most molecular bonds have the above heat-work 
transformation property, which models a chemical bond as a "'centrifugal 
heat engine"' . Expressions are derived for the standard Gibbs free energy, 
enthalpy, and entropy where the bond coordinates need not conform to a 
non degenerate Boltzmann state, since bond breakdown and formation are 
processes that have direction, whereas equilibrium distributions are derived 
when the Hamiltonian is of fixed form, which is not the case for chemical 
reactions using localized Hamiltonians. The empirically determined Gibbs 
free energy from a known molecular dynamics simulation of a dimer reaction 
2 A ^ A2, accords rather well with the theoretical estimate. A relation con- 
necting the rate of reaction with the equilibrium constant and other kinetic 
parameters is derived and could place the commonly observed linear relation- 
ship between the logarithms of the rate constant and equilibrium constant on 
a firmer theoretical footing. These relationships could include analogues of 
the Hammett correlations used extensively in physical organic chemistry, as 
well as others which are temperature dependent. One prediction of the prin- 
ciples developed here is that the equilibrium standard reaction free energy is 
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more dependent on the height of the intermolecular potential than its depth, 
so that the sign of the AG® can change for varying barrier height with fixed 
well depth, which may appear counter-intuitive. All the above developments 
can be tested directly in simulations and therefore provides a fertile ground 
for further research with significant implications on how standard states are 
determined in relation to the direction of chemical reaction. This work treats 
the molecular bond using standard thermodynamics as if it were a system, 
and it is anticipated that with the advent of single-molecule science and exper- 
iment, that might be one way in which molecular statistical thermodynamics 
would develop. 

1 INTRODUCTION,DEFINITIONS AND MODEL 

The Bohr correspondence principle can be used in the standard sense to de- 
rive the quantum description from the mainly classical description given here 
by appropriate mapping of the classical variables to operators and their inner 
products (norms). The standard states used is described below. Normally, 
the standard state is considered '"arbitrary"' and set to unity and the ac- 
tivity coefficients are relative to this chosen concentration state cq (e.g. unit 
concentration or pressure, such as moles per litre, Bar or Atmosphere); two 
primary concepts are implied in such an arrangement: (a) the assumption 
that the chemical potential expression ~ kT In j^i^ refers to setting to zero 
the work of the intermolecular forces (perfect fluid) relative to expansion from 
the reference state where the activity coefficient ji = 1 and where all 7's are 
dimensionless and only differences in work are detectable and (b) in other 
regimes, the isothermal work done is ~ kT In ^272) ^ or ^ e transition between 
two states at concentrations c\ and C2 respectively. The cqJq term is absorbed 
in the standard chemical potential (j,f where [ii = fJ,f+ kT In -^j^- The absorp- 
tion is interpreted as follows where the activity aj of a component is aj = Ciji 
for systems where the standard state is at zero density, leading to an apparent 
singularity. The chosen standard state (which we define and provide values 
for) used here nf(T) will be at infinite dilution or zero concentration (q — > 0). 
For the species (chemical) reaction Y2i=i u iA% = the standard free energy 
AG e "(T)is AG fr (T) = Ya=\ u if4( T ) + cannot be functions of density or 
any other variable other than temperature, where the conventional chemical 
potential for the work done from this zero concentration state to the one at 
d is 



w = M f+fcTln^i (1) 
dCi 

where <5q — > for the total isothermal work done for the standard state at 
Q = by simple integration of the work given by kT J Clearly eqn.© 

is not uniquely defined for the present and tends to a singularity. We can 
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rescale the potential and write 



where 5ciN = 1 for all <5cj.We choose 5ci to be arbitrarily small such that ji — > 
1 as q — ► 0. This limit accords with the Debye-Huckel theory of electrolytes, 
and in most other conventional descriptions. Then 



^{T)=^{T) + kT\n^^ ] (3) 
Setting the ^ of Q to be the same as ^ of (/Uj = 



tf+kTln^ = fif +kT In Ni + kT In 

dot ociNi 

or 

/i<(T) = /if + kT In iV< + A:T In ^ 
leading to 

£ , (T) = £(T) + T.Q i (4) 

where Qi = kTlnNi is an arbitrarily large number; the subscript i refers to 
the '"base"' material state for our development of reactions founded upon 
a basic unit (e.g. a proton or an elementary particle from which all other 
material states are referred to ). For instance, for the reaction raA ^ A n the 
standard state //^' (T) for 5ci — > must be of the form 

/# n (T) = (tiTQa„) + ni&(T) + 5\ n {T) (5) 

where 8^ (T) is the Gibbs Free energy of Formation relative to state i; ex- 
plicit expressions will be provided for the computation of 5\ (T) for n = 2 . 
Further define (nfi + 5% n (T)) = ^ (T) so that M * (T) = ^ (T) + nTQA n 
The above therefore is the form of the activity coefficients 7 used here. The 
other important presupposition used here and also implicitly in the literature 
that must be emphasized is that in writing down this form of the chemical 
potential, only problems about thermodynamical equilibrium are solved. It 
is an assumption to extend these forms to nonequilibrium regimes. From the 
Gibbs criterion, it follows that the standard state of the chemical potential 
for the potential [ij = /J%(T, c — ► 0) + ln^jC^ is such that 

M 

- AG fr = - t4 u J = kT ln K e ( 6 ) 

J'=l 
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so that nf(T,c -» 0) may be chosen so that © obtains, where K e always 
refers to the equilibrium constant. Since c — > is a limit and a value, the 
vant' Hoff equation also follows from such a limit due to cancellation of the 
specific entropies [U p. 182-183] leading to the equation 

d In K e _ AH® 

dT ~ ~kT>' W 
From the Gibbs criterion, we have 

M M 

- AG fr = - Wi = - E = ln K * ( 8 ) 

i i 

so that /it^' has a recurring term nTQ which cancels out in equilibrium ther- 
modynamics; we therefore by convention ignore this contribution by writing 
the chemical potential for substance i as 



^ = tf(T) + kT\n^ (9) 



From |H1 

-AG = = ^ vail = ^2 vnl, L = Y 



M M M 



which implies that — AG* = lnif e = — Vi/if. We note that 

M M M M 

^2 »n4* = E + 52(mTQ>i = E ( 10 ) 

i i i i 

because (niTQ)vi = 

Deduction 1 Subject to the truth of Gibbs' equilibrium criterion, the stan- 
dard chemical potential of species i fif is its free energy at zero density ( 
i.e. the free energy of a single particle entity coupled to a thermal well) and 
the chemical potential ^i(T, Cj) at any other state may be written /ij(T, Cj) = 
nf(T) + kT In cm . 

Deduction 2 From the standard state of any other composite species 
relative to the base species may be computed as //^ = rifi^ (T) + 5^ (T), 
taking into account the free energy change to that state. 

Deduction 3 Since the temperature parameter is specified, any species Ai has 
single species thermodynamics i.e. that of interactions with a thermal reser- 
voir at zero density, where the interactions are solely external forces between 
the reservoir and the particle, and the equilibrium thermodynamical proper- 
ties are the time average values for as long as the particle exists, and the only 
exchange of energy between the molecule and the external universe is via the 
heat reservoir. 
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Postulate 1 The Gibbs' postulate that the ensemble average equals the time 
average obtains, so that the species energetics are determined by the average 
of its motion in the phase space of associated with it. Moreover, since these 
are all single particle species that react or interact, they need not assume 
Gaussian-type distributions for non- degenerate states for the energy terms as- 
sociated with the reaction where there is an effective change of the Hamiltonian 
associated with the motions, as when a switching Hamiltonian is used. 

Definition 1 A molecular species A n (composed of n elementary A units) is 
an entity whereby the intermolecular or internal forces can be unambiguously 
distinguished from the external forces. 

Deduction 4 From Deduction^ if the single species is thermalized, the en- 
ergy exchange with the reservoir can only be the result of the external forces 
of the reservoir applied to the molecular species and vice-versa, and hence the 
energy distribution of the system would conform to the average energies of the 
various types of external composite motions (i.e. rotational and translational) 
relative to the composite body A n . 

Deduction 5 Since the species is thermalized at any instant, the energy equa- 
tions combine at all times the dynamics due to the internal and external forces, 
and the stochastic average of any "'static"' thermodynamical state variable is 
the result of the application of Postulate Reacting species exist for a finite 
lifetime during which the internal and external forces are distinguishable. 

From Deduction©, it is possible to couple stochastic and dynamical laws 
to determine the time evolution, e.g. conservation of angular momentum and 
energy equipartition. 

Postulate 2 Insofar as the single composite body (e.g. molecule, particle, 
etc.) is an 'external body' with respect to the thermostat, its thermal energy 
modes ( e.g. rotational and translational motiion about the centre- of -mass 
(CM)) are determined by the laws (quantum or classical) regulating these 
motions. 

Definition 2 A physical or chemical species species which conforms to the 
above description for standard states, Deductions midland Postulates (GK^ is 
thermomechanically coherent. 



2 DESCRIPTION OF A THERMOMECHANICALLY CO- 
HERENT REACTING PARTICLES AND MOLECULES 



For the reaction 



A n + A m <— > A; 



m+n 
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Figure 1: Diagrams depicting the rotational axes and '"bond"' formation j due to 
internal forces 

a general scheme is depicted in Fig. 1)2.1(1 with the appropriate coordinates. 
The simulation result reported here is when n = m = 1 for a simple dimer 
reaction. This model is extended by considering product formation about 
any bond j, depicted here along the z coordinate. IF masses A n and A m were 
permanently bonded along j, then the Hamiltonian for the system would 
permanently feature the intermolecular potential energy (p.e.) between sub- 
stances and the mean energy (for quantum systems) or equipartition of energy 
results from the density P having form 



where g represents the density in the total Hamiltonian space (p,q), which 
is not the same as the internal coordinates of the molecules, where if the spa- 
tial coordinates R = q in the system subsumed that about the bond, then the 
vibrational motion would exhibit the equipartition result, but for a reacting 
system with bonds that break (approximating a switch of potentials) , the 
phase density is a function of the total system density and such equipartition 
results would not obtain. For instance, Fig. (|2.1| ) shows the total internal en- 
ergy Ei n t distribution of the dimer in our simulation, where Ei nt = ^r 2 + V(r) 
where /i is the reduced mass, r the intermolecular distance, and V(r) is the in- 
termolecular potential (approximated as a Harmonic potential, where the end 
points which are curved in the actual potential was ignored). If there were no 
bond-beaks, or if they were very rare, then one would expect a Boltzmann-like 
distribution for the energy density, or the corresponding equilibrium densities 
given by QM, but this is not the case here since a real reaction is occurring 
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where the internal coordinates are not the immediate phase space coordinates 
of the entire system Hamiltonian TL; for very small values of the rate constant, 
the masses at either end of the bond j would be " 'quasi-ergodic" 'in that it 
would cover a fair portion of the phase-space of the internal coordinates, lead- 
ing to a density P given by the form (|12|) above. The internal energy is very 
approximately micro canonical in one arm of the distribution, TL{T) = E for 
the energy E and temperature T for the dimer studied here. This result should 
be contrasted with the many theories that presume that the reacting coordi- 
nate or internal coordinates have Gaussian or equilibrium energy distributions 
e.g. |H [3] or that the principle of local equilibrium obtains for these coordi- 
nates in the sense that Gaussian or equilibrium averaging is performed [3]. 
In the Eyring transition state theory, the vibrational coordinate connected to 
the '"bond"' does not exist;the current mechanical description accords with 
this insight. This result should be contrasted with advanced theories of bond- 
ing which all presume a Boltzmann or standard quantum probability factor 
of the form exp — TC(p, q)/kT (2HB1E] for various subprocesses in kinetics. In 
[21 p. 3318], the Maxwell-Boltzmann density is used to derive the polarization 
propagator in the binary collision approximation. The point here is that there 
is no reason to suppose this is necessarily true, since the relaxation times of 
other processes would determine whether or not such densities apply; for the 
total internal energy (vibrational energy) of the bond, Boltzmann-type den- 
sities might perhaps apply for relatively very slow reactions which have time 
to cover the phase space concerned for the potential of that configuration. 
The situation appears even more precarious in Aroeste's [31 p. 6 7] development 
where '"...we now further assume an equilibrium distribution in internal states 
of the react ants... " ' where in the simulation result here, the total internal en- 
ergy of the bond does not conform to these assumptions. In a recent review 
5, eq. (l),p.68] of standard Gibbs energy of association for large (protein) 
molecular fragments with ligands of form 



which is of exactly the same form as our dimer reaction, the standard reaction 
free energy is a complex function of Boltzmannized configurational integrals 
of form 



which includes the internal coordinates, where the complex '"contains six 
degrees of freedom that represent the residual translational and rotational 
motions of the bound ligand"' with little vibrational coupling, since it is as- 
sumed that "'the effect of the ligand's translational/rotational motions on 
either species' internal vibrational motions are very small. "'. The standard 
Gibbs' free energy is then given by expressions of the form 



A + B AB 



(13) 




(14) 



AB = - RThl {&* 



) 



( 



Zn,abZn,o 



Zn,aZn,b 



) 



+ P° (AV AB ) 



(15) 
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where R is the gas constant, T the temperature, C° the standard concentra- 
tion, the Z's are the various configurational integrals for solute solvent and 
dimer interactions, and P° {/SVab) is a pressure-work term. In the model pre- 
sented here, there is continuous coupling of the " 'external" ' modes with the 
internal intermolecular forces. There seems to be a strong possibility of elabo- 
rating the method presented here in rudimentary form to such systems where 
further generalizations (and perhaps some corrections due to oversight at this 
preliminary stage ) to what is given here would be required as more informa- 
tion and experimental facts are discovered. The science of thermomechanics 
and continuum mechanics [5] all use the standard heat and work terms of 
thermodynamics , together with various dissipation principles to couple pro- 
cesses and to describe the non-equilibrium properties. The use of tensors and 
other mathematical tools are in a very advanced stage so that all kinds of 
cross-effects connected to the geometry of the system (and symmetry) and 
coordinates of all the particles can be adequately accounted for. The work 
here focuses instead on the concepts of a new heat term and its connection to 
geometry, and doubtless the techniques currently available would be able to 
cast the concepts here into even more compact form. For instance, Lexcellent 
et. al [7j uses standard heat terms of equilibrium thermodynamical theory 
coupled with other terms within a continuum; others concentrate on an amal- 
gamation of basic linear thermodynamics description [SJ p. 2 10], where the 
laws of elasticity are derived from the Gibbs potential [El p.213,eqn.5]. Oth- 
ers are able to couple motions p.l84]where the actual interior forces and 
the motions and power dissipated can be described p.l79]but the difference 
between mechanical energy and thermal energy in relation to geometry and 
an external thermal source are not of immediate concern. 

2.1 Description of dimer and coordinates 

In the simplest case, the molecule can change "'shape'", here defined in terms 
of the moments of inertia about the 3 principle axis, with the change de- 
scribed by three orthogonal coordinates in the direction of the 3 Orthog- 
onal axes used to describe the molecule r x j,r y j,r Zt i t j, (i = l....r) where 
ri = [rx,li fy,Ii fz,i,l\- Obviously, one must choose the coordinates convenient 
for the system of interest in the general case. As in Fig.Q, the "'bond'" 
is along the z axis; if there is free rotation about the axis, then about 
bond j, i > 2 for r z i j and the moments of inertia Ij are functions of such 
changes of shape , i.e. for the simple case where there is no cross-coupling 
, we might write Ij = Ij(rj t i), (j = x,y,z; i for z coordinate only) . These 
moments of inertia pass through the CM of the product species, where if 
ly = J2 S m s Xi jS Xj tS , then = (i ^ j) where i, j are the coordinates e.g. 
X x , s = x s ,X y ^ s = y s ...etc. Concerning thermal energy, for quantum theory, 
the high temperature rotational partition function has the form |1U1 p. 78] 

/ 2 7 /-p\3/2 / j j j \ 1/2 
f rot where f rot = ( h z j ( - xx J y zz ) with the corresponding thermal 
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Figure 2: Non-Boltzmann distribution of the total internal energy of the dimer in 
this equilibrium simulation 

energy Ef given by 

i.e. for classical and quantum theories, there exists a function E t h errna i = 
Ety(T, I ,ri,TH) where TY refers to theory type (either classical or quan- 
tum) and TH always refers to the type of thermal energy (vibrational, trans- 
lational etc.); for translational energy about the CM, E tra ns = f° r each 
orthogonal direction; the same results obtain quantum mechanically. For 
classical rotation about the z-axis (free rotation) we would have an energy of 
{kT)/2 per independent rotation axis, leading to kT for the z-axis, for the 
free and independent rotation of the two reactant portions A n and A m about 
this axis. The situation is much more complicated for TY=quantum; for free 
rotation, the partition function is of the form ff ree = (I rec ;T) 1 ' /2 , I re ^ = 

I tnj™ j w h ere ry j s a constant, and the partition function for internal rota- 
tions with various activated states of angular rotation have the form /j. r . = 
X^o ex P — ( u i ~ u; o)/(kT) but the exact form of Matthieu's equation is not 
known for the solution of the Schrodinger equation for hindered rotations |10[ 
p. 89], and so several semi-empirical methods involving constants like 7 above 
have been devised by Pitzer jlll I12j and others where the total partition 
function for rotation ftot can be written ftot = ffreefi.r. yielding the thermal 
energy from (|16|) . The theory here is independent of TY; from the Bohr prin- 
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cipal, one can start from classical considerations and from there derive the 
quantum expressions by applying the associated operators. 

2.2 Determination of work-heat energy balance with molec- 
ular distortion 

For rotational motion, consider an element of thermal energy due to rotation 
about the given (z) axis given by the form 

E TY = E TY (T, I z , TH = rotn, z) (17) 

Then distortions to the molecule by the change in I z , in conjunction with 
the application of the external thermal reservoir would cause a change of 
the thermal energy by amount Ety(T,I z + 5I Z ,TH). On the other hand, 
prior to the application of the thermal reservoir interaction, the mechanical 
change of energy of rotation (with conservation of angular momentum when 
no external forces are acting on it) is E^cchilz + SI Z ,TH = rotn, z). This 
mechanical energy is derived directly from mechanics, (quantum or classical) 
without statistical mechanics since the body is a purely mechanical object at 
this interval of time. On the other hand, since the object must exist in a 
temperature field, and be in equilibrium with this thermal field and thermal 
reservoir, extra work of rotation <5u>TH=rotn,z m ust be provided externally to 
the body to make it compatible with the thermal environment which would 
have the rotational energy Ety(T, I z +SI z ,TH). Here, I z refers to the moment 
of inertia about the z axis (i.e. the dimer is rotating in a plane perpendicular 
to the z-axis) There is conservation of angular momentum during the dis- 
tortion, and the tangential impulses are applied by the reservoir to increase 
the kinetic energy (K.E.)to Ety(T,I z + SI Z , TH). The analogy of this ma- 
chine is the Carnot piston device, where if there is distortion of the cylinder, 
the Ety(TH = transl. and internal potential energy) thermal energy of the 
working substance is altered, and so the piston would have to be altered in 
position (by doing work) to maintain the same thermal energy density, where 
the work is 5w p = —PdV ; (here P is the pressure and V the volume) and 
thus we are lead to the following definition by the above analogy: 

Definition 3 An ideal centrifugal heat engine is defined as one where the 
work transfer increment 5w is the difference in energy due to the loss of ro- 
tational (mechanical ) energy, due to the change in shape or distortion deter- 
mined by the coordinates R z and the prevailing thermal energy in a temper- 
ature field at R z + dR Z; where there is conservation of angular momentum 
during the distortion prior to the application of coupling moments orthogonal 
to the axis to restore the rotational energy to that of the rotational thermal 
energy at R z + dR z . 

Likewise, we can define a rectilinear heat engine thus: 
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Definition 4 An ideal rectilinear heat engine is one where losses in thermal 
energy along an external boundary coordinate during the application of an 
internal force by the system where there is total conservation of linear mo- 
mentum can be compensated by the application of an external force in the 
same direction, and where the work done by the external force is the differ- 
ence in energy due to the loss of thermal energy along the coordinate and the 
prevailing thermal energy at that same temperature. 

For rotational motion, the mechanical, non-thermal rotational energy (denoted 
EMech,Rot,z{I(Rv), fi) would suffer loss 5EMech if SI Z > because of conser- 
vation of angular momentum |13| sec. 11. 6.1, p. 253] and the net work dVFthat 
would have to be done to enable energy balance for an ideal centrifugal heat 
engine is 



For instance, if only the z coordinates are involved then a simplified ex- 
pression is 



where T.E denotes a pathway along that of thermal equilibrium, , Fj, z is the 
resultant internal force (due in part to the interaction with the thermal reser- 
voir so it is not the pure inter particle force that the molecule would experience 
in the absence of thermal interaction) acting along the coordinate R z for the 
simple example where cross-coupling of forces has been eliminated. Eqn. (|19j) 
is a type of path differential, where the change - E g^ ch is along a thermal 
path that is consonant with equipartition of the thermal energy path at the 
specified temperature, as provided by quantum statistical mechanics. For 
simple systems, it can be shown by direct computation that the work transfer 
(heat absorption) is at a maximum if (|18|) is utilized along each stage of the 
path 5R Z . It is suggested that non-rotational motion interconversion such as 
what obtains above for rotating systems is the essential method of transfer 
for work- heat balances in normal thermodyanical systems, which leads us to 
the following: 

Conjecture 1 The rectilinear heat engine refers to the standard work-heat 
transformation of a fluid due to P — V ( Pressure- Volume ) changes, such as 
is routinely discussed in thermodynamical treatises. 

Theorem 1 The change of the Gibbs energy between two states AG due to 
the transition between two coordinate points R Z ,R Z can be written 




(18) 





(20) 
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Proof. There is an interconversion of energy between the external thermal 
angular momentum K.E. and the internal energy along R z brought about by 
an external force so that there is no net heat transfer (pure work of Gibbs 
energy conversion at constant temperature). Hence (|19|) yields a change of 
the Gibbs energy between two states and the above results • 

Hypothesis 1 The particle (including chemical) reactions found in nature 
are to a good approximation thermo-mechanically coherent if the reactions 
occur in an equilibrium temperature field. 



3 TEST OF HYPOTHESIS 



In the following subsections , the model of a reaction is introduced, which 
is applicable to all elementary second order reactions, and the general equa- 
tions, derived as theorems are applied to the above hypothesis to determine 
the standard energy states (standard Gibbs free energy, entropy and enthalpy 
of the reaction) . Estimates of the Gibbs energy is made based on the assumed 
shape of the molecular shape distribution function defined here; the determi- 
nation of the function constitutes a challenging but solvable problem which 
will be pursued in forthcoming investigations, and the current theory can be 
verified from the said function. The estimates give results which are already 
in good quantitative agreement with the Gibbs energy determined directly 
from the equilibrium constant. 



3.1 The model 

The dimeric particle reaction is 

2A - A 2 (21) 

(n, m = 1) above the supercritical regime of the L J fluid. The model resem- 
bles that of ref. |XI| except that a harmonic potential is coupled to the prod- 
ucts to form the bond of the dimer whenever the internuclear distance reaches 
the critical value rj between two free atoms A and no virtual molecules are 
formed by labeling or coloring the product. Hence the internal states of the 
molecule can be probed, from which the standard states may be determined 
according to the theory developed here. In the current study, the potentials 
as given in Fig. 01 are used, but other configurations are possible. In this 
model, the dimer A 2 is defined to exist the moment it passes the threshold at 
rj where the u mo i intermolecular potential is activated. The MD mechanism 
for bond formation and breakup is as follows. The free atoms A interact with 
all other particles (whether A or A 2 ) via a Lennard- Jones spline potential and 
this type of potential has been described in great detail elsewhere |15j . An 
atom at a distance r to another particle possesses a mutual potential energy 
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Figure 3: Potentials used for this work 



ulj where 

for r <r s 
for r s < r < r c 

and where r s = (26/7) e <j ^Hl- The molecular cut-off radius r c of the spline 
potential is such that r c = (67/48)r s . The sum of particle diameters is a and 
e is the potential depth for interactions of type A-A or A-A2. The constants 
ay and bij where given before ^5] as 

ay = -(24192/3211)e/r2 
6 4i = - (387072/61009) e/r 3 s 

The potentials for this system is illustrated in Fig. |31 Any two unbounded 
atoms interact with the above ulj potential up to distance rj with energy 
e = ULj(rf) when the potential is switched at the cross-over point to the 
molecular potential given by 



ulj = oy(r - r c ) 2 + bij(r - r c f 
ulj = for r > r c 



u mo l(r) = u vib {r)s(r) + u LJ [1 - s(r)] (22) 

for the interaction potential between the bonded particles constituting the 
molecule where u V ib(r) is the vibrational potential given by ea. (|2*3*|) below 
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and the switching function s(r) has the form given by ea.(|24[). At regions 
r < r sw,s(r) — » 1 implying u{r) ~ u v ib(r), i.e. the internal force field is es- 
sentially harmonic for the molecule and at distances r < r sw ,u(r) ~ u^j, so 
that the particle approaches that of the free LJ type. Concerning the mech- 
anism for the switching, in quantum mechanical kinetic descriptions, switch 
mechanisms are frequently used for describing potential crossovers|16|. but 
from a classical viewpoint one can suggest that the inductive LJ forces due 
to the particle potential field (with particles having a state characterized by 
state variables s^j) causes the internal variables at the critical distances and 
energies mentioned above to switch to state sjh when another force field is 
activated for the atoms of the dimer pair. State sm reverts again to state slj 
at distances r^. The following values of potential parameters were used here 
(Fig. El): 

uq = — 10, ro = 1.0, k ~ 2446 (exact value is determined by the other input 

parameters),?! = 100, r/- = 0.85, rj, = 1.20, and r sw = 1.11. 

The intramolecular vibrational potential u v n,{r) for a molecule is given by 

Uvib(r) =u + ^k(r - r ) 2 (23) 
The switching function s(r) is defined as 

s(r) = ^—^n (24) 



1+ 



where 



s(r) —5-1 if r < r 
s(r) — > for r > r 



sw 
sw 



. The switching function becomes effective when the distance between the 
atoms approach the value r sw (see Fig. © A molecule is formed when two 
colliding free particles have the potential energy u{rf) whenever r = rj at 
the value indicated above. This value can be defined as the isolated 2-body 
activation energy of the reaction. A molecule dissociates to two free atoms 
when the internuclear distance exceeds (which in this case is 1.20). Thus, 
the dimer A2 exists the moment it passes the threshold at when the u mo i 
potential is activated, and the distribution of the dimer state is along the 
segments a = (0 < r < rj) and b = (rj < r < r^) leading toc = aU6=(0< 
r < Tj,). 

LJ reduced units are used throughout this work unless stated otherwise by 
setting a and e to unity in the above potential description. The relationship 
between normal laboratory units, that of the MD cell and the LJ units have 
been extensively tabulated and discussed 15 and will not be repeated here . 
For the system simulated here with the potentials depicted in Fig. ©, the 
switching function is operative upto r^, the distance at which the molecule 
ceases to exist, and where the atoms which were part of the molecule interact 
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with the potential ulj like other free atoms;bonded atoms interact with other 
particles , whether bonded or free with the u^j potential. The point rj of 
formation corresponds to the intersection of the harmonic u v n,(r) and u^j 
curves , and their gradients are almost the same at this point; by the Third 
dynamical law, momentum is always conserved during the crossover despite 
finite changes in the gradient. Total energy is conserved since the curves 
cross, and errors can only be due to the finite time step per cycle in the 
Verlet leap frog algorithm, which would cause the atoms to be defined as 
molecules at distances r < r/. Similarly at the point of breakup, there is a 
very small (~ 10 -4 LJ units of energy) energy difference between the LJ and 
molecular potentials despite using the switching function in the vicinity of 
the region to smoothen and unify the curves; the small energy differences 
at the cross-over points are less than that due to the normal potential cut- 
off at distance r c where the normal (unsplined) LJ potential is used in MD 
simulations. A new algorithm which is general in scope was used to conserve 
momentum and energy in this study and the thermostatting was at the ends 
of the MD cell, as would be the case in most real physical systems. If E p (r) 
is the interparticle potential (energy) and E m (r) that for the molecule just 
after the crossover, the algorithm promotes the particles to a molecule and 
rescales the particle velocities of only the two atoms forming the bond from 
v i t° v 'i (i = 1; 2) where v'; = (1 + a)v; + (3 such that energy and momentum 
is conserved, yielding (3 = ~ a ^ T (^^^) V2 ^ (f° r momentum conservation) and 
energy conservation implies that a is determined from the quadratic equation 
a 2 qa + 2qaa - A = with a = (vi - v 2 ) 2 ,q = and A = (E p - 

E m ) where empirically, success in real solutions for a for each instance of 
molecular formation is 99.9 % and 100% for breakdown-where the A value in 
this instance is very small ( ~ 1,0 x 10" 4 ). This new algorithm coupled with 
shorter time step ensured excellent thermostatting. The y and z directions of 
the MD cell have length 1/16 each (cell units). This shape is chosen because 
it is intended that future simulations will concentrate on imposing thermal 
and flux gradients along the x— axis, which would allow for more accurate 
sampling of steady state properties about this axis, conversion units The 
shape of the potentials and switching mechanism used here is surprisingly 
similar to discussions of the charge neutralization reaction given in ^Hj 

K+ + r -> K + I (25) 

except that these discussions do not explicitly mention the crossing over of 
the KI and K + I _ potentials at short distances (high energy), although there 
is reason to suppose that such processes may well occur. It seems very feasible 
that reactions that have an electro-magnetic force law component in the tran- 
sition state may well exhibit some form of loop pathway such as the model 
used here. 



15 



3.2 Energy distributions 



Since this is an equilibrium system, Postulate Q leads to the dimer existing 
along a series of states about segment c above. Defining r n = n5r m where 
5r m is a fixed grid interval, we may determine the probability function by 
ensemble averaging at low particle densities (p* = 0.03 to 0.08) by binning 
the number of occurrences N Tn at each of the intervals centered at r n = nbr m 



where M = mt(rb/5r m ). In the limit, P r (r) exists by Postulate (|T|), As the 
dimer evolves about the r coordinate, work dw is done on the dimer until 
the single molecule has a characteristic equilibrium temperature T; all dimer 
properties therefor, including the heat or work transfer characteristics must be 
averaged by P r {r): by the Gibbs' postulate, the standard free energy AG fr (T) 
(determined at p* — > 0), must also conform to the general deduction below 
for fixed temperature: 5G = 5H — T5S and since P = 0, 5G = 5U — 5q. 
The total heat absorbed by the reservoir — 5q must be computed, as well 
as the total energy over the entire trajectory; both of these quantities are 
determinable from basic MD simulations. The thermodynamical path chosen 
on the single-particle micro system can be a work or heat exchange pathway, 
provided the same end-state is arrived at for both these paths since G is a 
state function. We choose a work path in our example, rather than compute 
the quantities above. The height of the potential curve at rf is 17.5153 and 
this is a work term. From Definition (^Q) and Deduction ©, we can view 
the particles and their aggregates as purely mechanical systems, thereby the 
'"bonds'" and other internal dynamics, without invoking any approximations 
whatsoever, the example of the dimer reaction assumes classical mechanics 
(TY=classical). The computations were over a range of temperatures T* = 4 
to 20 so that mean values of AG can be determined over this range (where 
T* = kT). The average value of K c = j^M ([X] is the number density of 
species X) is determined for densities in the range (p* = 0.03 to 0.08) and 
this average was taken as an estimate of K e , the equilibrium constant since 
7, the activity coefficient — > 1; the fluctuations are considerable within this 
range and no discernible reduction in K c value could be detected. At T* = 8, 
each mode would have from classical equipartition energy Tf, where Tf = 4 
(energy per mode at rf. From these specifics, we can deduce that the thermo 
mechanical coherent thermal energy- work exchange along a bond between 
two internuclear distances is given by Theorem Qbelow: 

Theorem 2 The thermal energy-work exchange AW^ between internuclear 
distances r\ and r2 is given by 



so that 



P r (r = r n ) 




(26) 




(27) 
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L' f 

Proof. For each of the 2 rotational motions, Tf = jfj, where L,I and T 
represents the angular momentum, moment of inertia and the total kinetic 
energy, with the subscripts indicating the position. When there is distortion 
of the r coordinate, then / varies, and the mechanical energy change in this 
case is dE M ech,Rot,z where 

^Mech.Rot.z = 2~^^ ^ ^ 

where R z = r the intermolecular coordinates. Since Tf = E ro t per mode, 
and 2ITf = L 2 , leading to 2IE ro t = L 2 so the increment of work dWc by 
substituting in ({25]) is given by dWc = —I~ x E rot dE and since AW^ = 
dWc, the theorem follows • 

Thus, there is a coupling of I and L because of the thermo-mechanical 
coherence property. By Deduction ((SJ, we can compare result (JJFJwith the 
thermalization requirement of Equipartition, where from 2ITj = L 2 we have 
2/4 = L 2 , and from (|19|) of the equilibrium thermal path with known tem- 
perature T, we have 

dW G = -AI' l dI (29) 

We note that E ro t is the equipartition energy or energy per mode (either the 
quantum or non-quantum energy, depending on which theory is considered 
more suitable). Averaging leads to thermodynamical properties, such as the 
one below. 

Theorem 3 The total change of the free energy about the bond trajectory 
AG mo i is given by 

= AW;*P(r 2 )dr 2 = AG mol . (30) 

Jr=0 



Proof. Since all dimer states are equally valid, the work done W r f on the 
single dimer (for all members of the ensemble involved in this motion) must 
be (from l|27|) and (J2EJ) ) (per mode) the time average or ensemble average 
value leading to the result above • 

The above AG mo i is the increment of free energy associated with the bond 
in its trajectory. In (|3U|) . for ant transition 5G = 5H — T5S for any transition, 
so that if 5q is the input of rotational energy caused by internal stresses, 
then the internal work done is 5W In , and 5q is the heat absorbed where 
Sq = a and 5H = 5W + Sq, then —T5S = —Sq, so that 5G = a + a — a = 
a = 5W. Alternatively, instead of expansion, we can compress the single- 
molecular system, and (|30JI is derived if we note that if the bond axis about 
z is compressed by doing work 5w, so that the two other variations S in the 
moments of inertia about the x and y principle axis obeys 5I XtV < 0, then the 
increased kinetic energy of rotation would have to be removed by amount —5q 
by the reservoir, where 5H = —5q+5w and since —T5S = T5q, then SG = 5w, 
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and integrating this leads to AW™ = JJ 2 Sw that can be inserted into ()3U[) . 
Finally some approximations may be derived. (Some exact expressions will be 
presented later. If the work done is approximately equal to the rotational 
energy dissipation, then SU rs Sw + 5q = or 5q » 5w, or T^S w — 5w or by 
integration TAS™ ~ AW™ and so 

T (AS) = T /" ' AS™P{r 2 )dr 2 » W^J (31) 



r=0 



or 



(AS) « ( ^ ] (32) 



r 

At present, the form of P is not available, since the simulation program de- 
veloped did not anticipate the need for this distribution, but it will be included 
in future investigations for the years ahead, and would provide key evidence 
for or against Hypothesis Qand the current formulation. But estimates of 
P(r 2 ) can be made for extreme cases in the current study where the real value 
should lie in-between the extremes; this is what is found; further, the extreme 
cases yield at the least semiquantitative results, meaning that quantitative 
results and corroboration of Hypothesis Q is expected from a detailed study 
where P(r 2 ) is known. It will be noticed that the total bond energy E to t 
defined as E to t = U m (r) + ^r 2 is very approximately constant (Etot ~ c) ,and 
roughly simple Harmonic in nature, with the defining equation being 

I* = Fint (Vint = • (33) 

This motion would be moderated with the actual shape of the potential, 
and the collisions caused by the heat bath in the radial direction; such dy- 
namics leads to the curve in Fig. ()2. 1|) . It is conjectured that the rate of the 
thermostatting might have an effect on the shape of the curve, as it would lead 
to the boundary conditions of (|33|) to be reset each time there is a collision. 
The several scenarios which can occur are listed below: 

(1) the molecule is formed at rj and then traverses to in such a way 
that P{r2) = constant and P(r < rt) = 0. The motion approximates 
(|33|) except there is no accumulation of phase density at the turning 
points of the curve, so all points have equal weight (consistent with a 
transition stage with no multiple reflection back to the vicinity of rj 
with Etot = c ). The array of conditions may be written 

i(n) = i{r f ) = i Tf , 

I(r 2 ) = I(r b ) = I n , 
"(I) _ M„^7 



W rf ^ = -2(E rot )\n{jL) (34) 
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The parameters are r } = 0.85, r b = 1.20,1/ = 0.36125, I b = 0.720. For 
the 2 rotational modes at T* = 8, E rot = 4.0 for classical systems, so 
that (El yields W^ (1) = 0.6862T* 

(2) since AG fr (T* = 8) from the simulations, we can adjust Ifa) (I ri = Irt 
is fixed) so that there is exact coincidence at AG fr (T* = 8); this occurs 
when rf = 1.1539 (instead of 1.20) and /(r&) = 0.6656 where r/ and If 

(2) 

are as in (1) above. The W rf function, denoted in this instance W rf 
is from (JUJ) W^ {2) = 0.61110T* 

(3) here we choose I ra t = — Inf-^j = 0.50, where = 0.60653 and 

(3) 

the W r . function in this case is 

3.3 Derivation of standard thermodynamic functions for dimer 
reaction 

For the following transition scheme, 



2A — > A2 (threshold at rf = 0.85) — > A2 (dimer phase space) 
State 1 State 1* State 2 (35) 

the molecule '"exists"' in the space (State 1* U State 2). 

Notation: P : Pressure; U, E : Energy function ; H : Enthalpy ; T,T* : 
Temperature; asterisks imply reduced forms, e.g. e* := (kT) = T*; 5X 
: variation of X; w denotes a work term, and q always refer to heat variables; 
subscripts denote whether the dimer-molecule (m) or atomic species (a) is 
being referred to; L refers to the latent energy of the particles (which usually 
cancels out in the expressions). 

An elementary treatment of the statistics of a '"bond"' that can be dy- 
namically broken and formed is given in Appendix (A.) 
State 1 — > 1* transition 
Since P -> 0, H = U. So, 

H a = U partides = Qe* x 2 + 2E L = ^e* + 2E L = Swi + % + 2E^j (36) 
where 5w\ = (no work). For molecules, we get 

H m = V 'molecule = <5W m + 2El + <?1* 

= ~ e * + E L + 6W m (37) 

where q\* = 2(^-) + f e * f° r the translational and rotational energy of the 
composite structure. El is the latent energy of the isolated atom at the 
indicated temperature. It is crucial to realize that SW m is the work done to 
'"fuse" ' the particles at thermal equilibrium into a single body (from the point 
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of view of mechanics previously described) and that after that process, the 
external reservoir only works upon a single external body that distorts itself, 
and the forces of interaction are through external forces; (|37j) represents the 
rotational and translational K.E. about the CM and the residual energy which 
is carried over from the constituent reactant species and the work e)Ty m that 
is required to fuse the reactants, where SW m = 17.5153 and this comes from 
the p.e. at rj given in Fig. ©. The change in enthalpy is therefore given by 

AHt^ =H m -H a = 5W m -j (38) 
and the entropy change AS is such that 

TAS^* = - \e* = -\e* (39) 
leading to 

AG^p = 17.5153 - ^e* - TAS = 17.5153 . (40) 

A heuristic interpretation of the " 'bond" ' free energy is possible. The z axis 
translational k.e. is 2.(e*/2) and when the bond is formed, one can still 
refer to the heat content of the bond as qtcmd = e* • Since 5H = £, max + e* and 
-T5S = -e*, the bond free energy change is AG* bond = (£ m ax+e*)-e* = Cmax- 
State 1* — > 2 transition 

Thus, for thermo-mechanical coherent systems, as defined in (|30[) . we have 

AG fr (T) = AG(T) = A6?i_i. + £ (j^ dG^j P{r 2 )dr 2 
or 



AG (T) = AGi-^i* + W r f (41) 



with W r f defined as in (|3Uf> . In general, however, 

AG^(T) = AGW + ^ 2 (dH - TdS)^j P{r 2 )dr 2 (42) 

. Ean (j41Jl will be considered in detail in the next section. We also define the 
following energy term for the two cases that we may encounter: 



Wrf' 



W r f from (J3(J|) if system is thermomechanically coherent, 

otherwise we set . (43) 

J? (dH - TdS)) P(r 2 )dr 2 
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Figure 4: Linear Plot of \nK e vs 1/T* 



3.4 Estimates of AG^(T) 

From (|41j) . and 1)34(1 .we derive 



AG®^(T) 
AG® (2) (T) 
AG fr(3) (T) 



AGi_ 
AGi_ 
AGi_ 



+ W, 



rf 



"(1) 



+ w rf 



<2) 



+ w, 



r/ 



"(3) 



17.5153 + 0.6862T* 
17.5153 + 0.61110T* 
17.5153 + 0.500T* 



(44) 
(45) 
(46) 



The results from the simulation for AS* 6 ", AH® and AG 6 " follows. Here, we 
can only discuss estimates for AG* based on the theory; the derivation for the 
enthalpy and free energy of reaction will be presented but no detailed discus- 
sion is useful at this stage until the P{r2) function is determined numerically 
so that a comparison between the theoretical predictions and the numerical 
values derived from simulation is made; this is a long term challenge. The 
van't Hoff equation (at constant pressure, P — > 0) is 



dhxK e __ A r H*(T) 
d(l/T) " R 



(47) 



A plot of lnir e vs 1/T* is given in Fig. @ where the plot is rather linear over 
the range of experimental temperature points; therefore the gradient would 
yield a fairly constant average value applicable for the temperature range for 
the standard enthalpy AH® . The overbar refers to these mean values over the 
temperature range (T* = 5 — 20). Since AG® is known, the average standard 
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entropy may be determined from N readings as 



N 



^^E A(y(T t Ag J <48> 

The theoretically less correct method is to do a linear fit to the equation 
(assuming that both the standard entropy and enthalpy is reasonably constant 
without prior justification) 

AG fr = AH e -TAS®. (49) 
The van't Hoff (v.H.) results using ()47|) are 



AHf H = 18.8 ±0.6 (50) 

AS% H = -0.47 ±0.06 (51) 

AG%, = AHf^-TAS%. (52) 

At T* = 8, the experimentally determined AG 9 " (from extrapolation to zero 
density) is 

AG fr = 22.4 ± 0.3 (53) 

The v.H equation determination (|52j) for the same temperature is 

AG% = 22.5 ± 0.6 (54) 

which is close to the experimentally determined value. The linear fit (/) 
estimate (|49|) for the entropy and enthalpy are 



AHf = 19.2 ±0.9 (55) 
ASf = -0.43 ±0.07 (56) 

Surprisingly, perhaps, (|55(l56|) are reasonably close to the v.H. results. Fur- 
thermore, from estimate ()32|) . we observe that ASf is close in value to the 
range (0.500 -0.6862) derived from (I44I46|) . 



3.5 Interpretation of results 

The values (AG 9 ^ 1 ^) from (|44B46|) plotted against the experimental curve 
given in Fig. © is essentially quantitative even for the approximations made 
here. AG 9 ^ 1 ) (Case (1)) represents a distribution of "'dimer states'" uniformly 
distributed, and is quantitative for 4 < T* < 10 (lower temperatures). The 
cases referred to in the superscripts have the following moments of inertia 
ratios 

= 0.502; (IfA = 0.543; (-f±) = 0.607 . (57) 
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Relative to the loop mechanism, for fixed I rf , the results clearly show that 
there exists an effective value of the moment of inertia determined by P(r 2 ) 
designated I rb where I rj> < I Tbt x (the arabic numerals refer to the Cases 
(1-3) )for optimal results if the theory is correct, and that this moment of 
inertia would change with temperature; there is a non-uniform distribution of 
steady state distances with an accumulation of distance density at distances 
r where r < r^; i.e. there is a preponderance of states to the left of the region 
of dimer breakdown. Indeed, at high temperatures, the probability of states 
x with internuclear distances r x — > rt would increase in the vicinity 5 about 
rj for arbitrary 5 so that I Tb < I rb) i leading to the observations. However, 
the overall changes are not very dramatic. In particular, once the effective 
I rf /I rb ratio has been selected to coincide with an experimental point (Case 
2 at T* = 8.0) by adjustment of the parameter to fit with the actual P(r 2 ) 
distribution function according to the prescription below (where the notation 
for W is as in (|30JI ) where 

- T * In r-f] = f b AW?P{r 2 )dr 2 (58) 

\hj fit J r =0 1 

then there is quantitative agreement about the whole range for T < 15. We 
have also to derive from eqns. (|5QI51|) the AH® and AS^terms formally for 
this theory, which must follow from the internal dynamics of the molecule. 
For theorems 0] and EJ the term AW^ can refer to either that of ()41|) OR to 
the integrated Gibbs energy up to coordinate r 2 , that is 

{AW£j from ()27|) if system is thermomechanically coherent, 
otherwise we set 
(dH - TdS) 



Theorem 4 The standard entropy is given by 

AS*(T) = - j2 ^ (AW^P(r 2 , T)) dr 2 (60) 

Proof. The fundamental thermodynamical relationship 17 :i p. 167, Sec. 7.2] 
AS® = -^r- \p, a may be applied to (jUj) written 

AG®(T) = Uax +Wf= I* AW^dr 2 P(r 2 ,T) + Uax (61) 

Jr=0 

and differentiating this equation with respect to T leads to the resulting the- 
orem • 

Equation (|60jl is testable from the point of view of simulations, so that an 
independent and explicit determination of AS® (T) is available, which can be 
correlated with the entropy derived directly from curve fitting to the AG® (T) 
as detailed above in the previous sections. We note the elementary fact that 
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Figure 5: Plots of the estimation of the standard Gibbs' energy for the three cases 
discussed in the text together with the simulation results. 



entropic values are intimately connected to probabilities, and so ()6U|) is another 
example; indeed if TdS = Yli=o e i^ n i where ej are energy levels p. 642] 
and drii are the changes in population then dS ~ YliLo( e id n i) /T which is 
of the same form approximately as ()6Uj) . As with theorem (^Q) the standard 
enthalpy can also be derived. 

Theorem 5 The standard enthalpy of reaction AH (T) is given by 

AH*(T) = i max + jP (AW;; |p(r 2 ,T) - T 9P{ ^ T) 

dAW^ \ 

- P(r 2 ,T)T—^jdr 2 (62) 
Proof. Either by subtraction (AH®(T) = AG® (T) + TAS e ~(T)) or the ap- 



plication of the Gibbs-Helmholtz equation of the form^ y ^ j = — ^ 

we derive from 1)611 160)) the above result • 

Again, (|62J) is directly testable if the P function were known. The above 
equation also explains why £, m ax ^ AH® because there are other terms in 
(|62[) that contributes as well. 
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4 FREE ENERGY-LIKE RELATIONSHIPS IN KINETICS 



This section includes a discussion of the Arrhenius rate law and its relation 
to the equilibrium constant. If we admit that simple collision theory (SCT) 
or the quantum version due to Eyring leads to Arrhenius type equations for 
elementary reactions which determines the rate law then we can state the 
following: 

Theorem 6 There exists a relationship between the rate constant k r ^ and the 
equilibrium constant K e given by 

lnk rA = B(T) +h\K e (63) 
where the B{T) term has the form 



B(T) = In 



4(T)exp 



W rf ' 



(64) 



kT 

Proof. Eqn. (|6*Tj) savs that 

AG*(T) =(, max +Wf ' (65) 
On the other hand, SCT writes the rate v as 

v = k{e r )n\n2 with k(e r ) = cr(e r )v r . (66) 
The above expression refers to the rate for the internal kinetic energy e r = 



^rr. where the n's are the densities of the reactants. Integration of k(e r ) leads 
to the form of the forward (1) rate constant given by k r ^\{T) p. 99] with 
the rate expressions 



v = k Ty i{T)n\n2 with 
k r i(T) = av rei (fi, T) exp - 



kT 



= AiCTjexp-^J (67) 

where an Arrhenius form is observed in the last line of (|67|) : a is the 
maximum impact parameter, v re i(/i,T) is the relative velocity of the particles 
and e* is by definition the threshold energy parameter, whereby all reacting 
molecules must have this energy along the line of centers for reaction to occur. 
The situation is not altered in Eyring's TST, where e* goes one step beyond 
it, in that e* represents the zero-point energy difference between reactants and 
activated complex |181 p. 103]. Hence, in this model too, e* is the threshold 
parameter. For TST theory, the form of the rate constant is 

^ (T) = f^ exp( -^ ) = ' 42(T)exp( -^ ) ' 
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The q's are the partition functions of the intermediate ^ and reactants A, B 
[Tsj eqn. 3.17,p.l05]. Both theories are of form ^4j(T)exp — (Jy). For com- 
pleteness, we shall augment the forms k r ^ (which does not alter the form of 
the standard expressions, and which may be neglected in most studies not 
concerned with the energies considered here) by a work term c due to the mu- 
tual intermolecular forces acting on the reactants, which is a '"real"' term 
since these forces exist, and which has the form c = c(nA,nB,Q) where 
n At n B are reactant densities, and fi refers to the physical variables (non- 
thermodynamical) such as the dielectric constant or the screening parameters 
. We are focusing on an elementary bi-molecular reactions here as an ex- 
ample. We can also apply a steric factor P so that the rate can be written 
kr,i = P(T)MT) exp — ( 6 fcy C ) ■ As nA,riB — > 0, we expect c — > 0. Write 
exp -(^) = C"(T), so that 

k r>i = P{T)C'{T)A t {T) exp - ( . (69) 



kT y 

Define A'^T) = P (T)C' {T) A^T) so that 

fc r)i =4(r)exp-^V (70) 

We interpret the pre-exponential factor Ai (T) to be that which obtains (when- 
ever P is known and specified exactly) whenever the particle density p — > 0. 
Identifying the e* threshold with £ max , the threshold energy for coalescence 
leads to 



AG (T) - W r f ' = i m ax = £*• (71) 

Since the equilibrium constant K e is given by — AG*(T) = kT dn K e , we are 
lead to 



kr,i = A-(T)exp 



1 -{AG*-W7f'} 



kT 



4(T)exp-^j.K e . (72) 

Taking logarithms on both sides of Q72|) leads to the result •. 

Theorem© is cast in the standard language of the " 'free energy relation- 
ship"' linearly relating for instance the logarithm of the rate constant with 
the logarithm of the equilibrium constant |17[ p. 961, Sec 27.5] where it is 
opined that the activation energy is linearly proportional to the standard free 
energy change for a series of compounds, leading to correlation analysis used 
in physical organic studies, fashioned after the Hammett and Taft equations 
which are essentially an empirically constructed set of fitting parameters [191 
esp. Chap 3]. There are two different ways of approaching Theorem Q and 
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1) From a physical point of view, of little interest to physical organic 
kineticists, (|63|) may be viewed as an extention of linear Arrhenius plots, 
where lnk r ^ is plotted against 1/T* to determine the activation energy where 
A[(T) is treated as a constant over the temperature range of the plot whose 
slope yields the activation energy. In this sense, the linearity of the log-log 
plots using is dependent on the stability of the B(T) function of over 
a large temperature range. Some comments are in order here; that lnA^(T) is 
experimentally rather invariant is a fact observed in many chemical kinetics 
studies, where the Activation energies are determined routinely by plotting the 
rate constant against the reciprocal temperature; in these estimates, lnyl-(T) 
is presumed invariant and accurate values are determined. For the cases 
considered here, W r f ' are linear functions of the temperature, so that division 
by the temperature factor would lead to a constant expression; actually it 
would be weakly dependent on temperature due to the P(r2,T) function, 
which is expected to have a slight temperature dependence. Hence, for the 
same reaction conducted at different temperatures, we can predict, subjected 
to the above assumptions, a somewhat linear relationship between lnfc ri i(T) 
and ln-ff e . This prediction can be tested in simulations as well, for the general 
applicability of (JHSJ- 

2) At any one fixed temperature, if the steric P and C'(T) function repre- 
senting intermolecular work is reasonably constant, then for any two similar 
reactants A and B, subject to the above conditions, we would expect some 
forms of linear log-log relationships when the " 'bond" ' structure is fairly sim- 
ilar (e.g. similar activation energy and reaction pathway in the transition 
state). The correlations would be found if various forms of linear relations 
are found for W r f ' for the same reaction with different substituent groups. 
In view of the fact that there is nothing unique that is specified by Theorem 
© concerning linear relationships, a great many types of correlations may 
be expected by the expansion of other functions that can be related to W r f ', 
such as ^max an d the masses of the particles constituting the molecule, since 
these affect the moment of inertia I , and the translational motion. Such an 
undertaking is a deep study in itself which will not be attempted here, suffice 
to say that Theorem © can provide a foundation for this undertaking. 

5 CONCLUSION 

Particle reactions moderated or subjected to a temperature field may be 
treated as a three dimensional (3D) mechanical body which interchanges its 
internal energy where these interchanges are due to the interaction of internal 
forces within the molecule and the thermal reservoir due to change of shape 
of the composite particle system or "'molecule'". The analysis above allows 
for the treatment of single-molecular thermodynamical systems. Advances in 
technology very recently are gradually bringing about the realization of these 
entities |2U1 1^ . It is possible to predict the standard states of such single- 
molecular species by measuring its internal motion at equilibrium where the 
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P[r2) probability distribution function serves as the key in the derivations and 
this probability function can be determined in principle from single-molecular 
studies. The linear "'Free Energy'" relationships are given a firm theoretical 
basis because of the derived relationship between the equilibrium constant 
and the rate constant for the elementary reaction. A deeper study would 
derive relationships similar in form to that provided by Hammett and Taft. 
It is predicted that the height rather than the depth of the potential well 
would determine the sign of the standard enthalpy and free energy.All these 
assertions can be tested, and future investigations would be directed toward 
this end. 

Appendix 

A Time Average Probability Densities From Threshold En- 
ergy Values 

The definition of the activation energy is that for r = v r = = -E^m (where 
Ekin = \^v^.)xi max is the minimum energy required to form a "'bond'". Dur- 
ing the formation of the bond at the threshold rt, the values of the internal 
kinetic energy has the range < Eki n , +00. We choose the time increment 8U 
to be small enough in the grid interval so that the following obtains: 

\5r\ = \ri + i — rj| = (3, where |£ijiv| < P (73) 

For E^in = a at rf when a molecule is formed, the work change bW a ^ within 
two adjacent grid sites is 



5W a>i = {W rf {r l+l ) - Wrfinj) + (74) 

at most by property (|73|): if the bond coordinate difference is < /3,then SW a ,i = 
0. We can envisage the intermolecular coordinate r% hopping by one increment 
(i.e. from i to i + i or i — 1). Define an interval SE in energy where E a — ^ < 
E a < E a + $y The probability density P r ,a( r i) f° r a series of individual trials 
j for the formation of the molecule to its disintegration about energy E a of 
interval 5E about the grid coordinate 8r is 

PraM = lim N' — ► 00 *• ^=i;;^ } with ( rN rb (75) 
Since a varies from to +00 we can define the total probability as 



P(r, T) = lim 



M 
1 

"fir 



OO 
OO 



E 



a=M&a 
a=0 



Pr.aiTi) 



^a=M5a 



=0 



PrjOtiVi) 



(76) 
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and (fTS)) must be equal to l|25|l if the Gibbs' postulate ((J) is correct for a 
sufficiently dilute system. The shape of the P(rz,T) is determined by the E a 
distribution, together with the internal forces of the bond (molecule) and the 
mean relaxation time of interaction of the thermal reservoir and the molecule. 
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